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Summary. Abstract We will discuss some specific applications to the rotation state 
and the shapes of moderately large asteroids, and techniques of observations putting 
some emphasis on the HST/FGS instrument. 



1 Introduction 

Although their name suggest their are point-like, asteroids are from long time well 
known to show variations in their lightcurves with shape and rotation [1] . Observed 
lightcurves can for instance be explained by spinning tri-axial ellipsoids, but even 
better by convex shapes in uniform rotation [2, 3]. Also the rotation period of these 
bodies seem to have some connection with their size. For instance asteroids larger 
than approximately 0.15km do not spin faster than ~ lOcycles/day (see Fig. 4). A. 
Harris [4] has suggested that this limit is not the tail of some statistical distribu- 
tion but does correspond to the limit of disruption of a gravitationally bound and 
cohesionless rubble pile. We will discuss in the following on the inversion of aster- 
oids lightcurves, on rotational state of asteroids in general, on possible figures of 
equilibrium for rubble pile asteroids, and illustrate how we can derive information 
on asteroids shape and size from high resolution observation with the HST/FGS 
interferometer. 



2 Lightcurves 



One of the first and most important source of knowledge on the asteroids physical 
properties was obtained from photometric lightcurve observations. These reveal the 
asteroid's rotation period, their brightness variation due to their non spherical shape, 
albedo spots, light scattering of their surface, ... Two typical examples of asteroids 
lightcurve are given in Fig. 1. In particular one can see from Sylvia's composite 
lightcurve that the brightness variation is periodic, and that there are two - almost 
identical - maxima and minima. Thus Sylvia's brightness variation can be well be 
approached by a tri-axial ellipsoid in rotation. So that analysis of several lightcurves 
obtained at different apparitions provide the pole orientation and the axis ratio of 
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the ellipsoid shape model. In the more general case however, lightcurves are not 
always so smoothly sinusoidal, but more irregular and asymmetric with additional 
extrema, ... (see the lightcurve of Hebe in Fig. 1 for an illustration). Promising 
results in the lightcurve inversion problem have been obtained recently [5]. 
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Fig. 1. Composite lightcurves of two asteroids from the Asteroid Photometric Cat- 
alogue [2]. Upper panel: (87) Sylvia observed on 1987 Feb. 3.3 {filled circles), and 
Feb. 6.3 {filled squares), by Weidenschilling et al. (1990). Lower panel: (6) Hebe 
observed on 1987 Jun. 18 {filled squares), Jun. 23 {filled circles), and Jun. 27 {open 
squares) by Hutton & Blain (1988). 



3 Rotation 

Considering a freely rotating rigid body (Euler's spinning top), integration of the 
Euler equations yields the orientation in the frame of the body of the instantaneous 
velocity fi, and the orientation - via the Euler angles (0, 0, tp) - of the body in a 
inertial frame [6]. As noted in [7] such a dynamical system has two - generally non 
commensurable - frequencies, so that the body never shows the same aspect in time. 
Although this exists in the solar system, there are however only a very few small 
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bodies that are known to clearly show (or suspected to show) such a spin state [8]: 
among them one comet P/Halley [9], and one asteroid asteroid (4179) Toutatis [10]. 
Moreover all these objects have relatively long rotation (spin) period. We shall see 
later the reason for this lack of (fast) precessing bodies. 
Let us remind that integration of the Euler equations: 

II Ull + {Is — h) 0^2 W3 = 
h d)2 + {h — Is) wi = 

+ (/2 — /i)a;ia;2 = (1) 

and 

• uji sill (,,' + ^'2 cos 'Ip 
sin^ 

; . • cos d (uJi sin lb +u)2 cos tb) 

tp = uJs — cosd <p = uJz ^-7 

sm^ 

= u)\ cos %[) — u)2 sin ip (2) 

gives the orientation of the asteroid principal axes with respect to an inertial frame 
at any epoch t [6, 7]. In a previous chapter by Tokieda we have also seen that 
the conservation of kinetic energy and angular momentum provides two integral 
relations: 

2E = Iiuii + 121^2 + I3 UI3 

= ll ljI + 7| w| + /| w| (3) 

The (xi,X2,xs) body-frame (see Fig. 2) is a right-handed frame associated to 
[xs,x,,xi). and with this choice of indexing we have put Ii > I2 > Is,- In gen- 
eral the Euler angles are given with respect to a inertial frame where the 2;— axis is 
- for commodity - aligned with the angular moment, and are understood as rota- 
tion, precession and nutation angles. In that case both angles 6 and ip are periodic 
functions of commensurable period. Putting: 

2 ^ (J2-Jl)(2gJ3-M^) 

(73-/2)(M2-2£;7i) 

and making use of the elliptic integral of the first kind if = J^^^i^—k^ sin^ u)~^^^ du, 
the period of the rotation angle ip is [7] : 



^-=^^^/ (/3-72)(M^'l2E/0 

while the period of nutation 6 is Tnut = Trot/2. Last, Landau & Lifchitz [7] have 
shown that the angle (p can be obtained as a sum of two periodic functions (p{t) = 
(pi{t) + 4>2/{t) where the period of </>i(t) is exactly Trot/2 and the period of </>2(i) is 
T', which in the general case is not commensurable to Trot. The latter period can be 
obtained from 2 tt Trot / /q^'°' 4'{i) dt. For celestial bodies one may prefer to express 



4 



D. Hestroffer and P. Tanga 



the Euler angles in the frame of the echptic J2000. In that case the two fundamental 
frequencies are mixed, and the Euler angles can be expressed as a sum of periodic 
functions of non-commensurable periods. Interestingly, none of the Euler angles is 
varying uniformly with time in the general case, which gave the name 'tumbling' to 
the asteroid Toutatis spin state^. In the case of Toutatis, two of the principal axes 
have approximately the same moment of inertia (72 ~ 13), so that the nutation 
is small {6 ~ 0) and the precession and rotation are circulating with non uniform 
velocity as shown in Fig. 3. Note also that, considering the readily inequalities 
2EI3 < < 2EI\, we can distinguish long-axis mode rotation (LAM, which 
is the case for Toutatis'* when Af^ > 2E I2, from short-axis mode rotation (SAM) 
{M^ < 275/2). This denomination reflects the fact that the instantaneous rotation 
axis is closer to the long axis, or the short axis of the body (excluding "medium-axis 
rotation" that anyway are not stable). 




Fig. 2. Euler angles {(f), 9, tj)) of the body frame {xi,X2,X3) given in an inertial 
reference frame {x,y,z). 

One can thus expect that rigid bodies that suffer collisions and /or external 
perturbation can be in complex rotational state (misaligned rotation). In the general 

^ See an animation on URL 
http : //www. star .ucl . ac .uk/ apod/ solar sys/ cap/ast /tout spin. htm 

* This is the reason of our choice of ordering Ii > I2 > I3, when the most often 
used one is to put I3 > I2 > Ii- With the presently adopted convention we only 
ensure that the rotation is associated to the Euler angle and axis X3 etc. 
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Fig. 3. Euler angles for (4179) Toutatis. Note the non uniformity of the circulation 
for angles and ip. 



case however the asteroids do not show such large free precession/nutation tumble 
or wobble over time span of typically a few days^. In contrast, it is observed that 
they arc spinning around an axis that 1) approximately coincides with the direction 
of maximal inertia, and 2) which direction is approximately fixed in (inertial) space. 
Asteroids, as the majority of celestial objects, are not infinitely rigid bodies but 
are deformed under stresses. When not in principal axis spin-state, there is a loss in 
kinetic energy: during each wobble period a portion of the stored stress-strain energy 
is dissipated in the asteroid's interior. Since the angular momentum is conserved and 
2T = / 1 is decreasing, this dissipation is balanced by a spin state that evolves 
asymptotically in time toward a rotation along the axis of largest inertia (/i with the 
notation adopted here). The timescaJe of this damping process has been analyzed 
by [11]. It depends on the spin rate f2, the shape and size D of the body, its 
density p, and of course on its rigidity /u: 

Tir.^lQ/{pKl{D/2fn^) (5) 

where Q is a quality factor expressing the ratio per wobble period of the energy 
contained in the oscillation to the energy lost. If the body has no rigidity it will 

* Precession due to planetary perturbation is not discussed here. 
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Fig. 4. Spin periods versus size of 987 known asteroids (from P. Pravec.) and 
damping time scales {dashed lines). A limit in rotation rate for large bodies has 
been drawn by hand. 



instantaneously align its spin axis, and if it is infinitely rigid there is no evolution. 
Considering values typical of asteroids, [12] alternatively gives the damping timescale 
by: 

Td ~ pV(17^-C>') [inlO« yr] 
where P is the rotation period given in hour, and D is the diameter in kilometer. 
Asteroids densities are (roughly) in the range 1-3, their shape gives Ki ~ 0.01-0.1. 
The value of the product nQ is not known and is still in debate. Burns & Safranov 
used a value of ~ 3 x 10* (cgs units), on another hand Harris [12] - by considering 
rubble-piles instead of non-porous solid rock and based on available estimates from 
an analysis of Phobos orbit - adopted a value of 5 x 10^^ Pa, while Efroimsky [13] 
states that the damping time would be 2 to 6 orders of magnitude shorter. Damping 
timescales with values adopted from Harris are plotted in Fig.4. One notes that 
either small asteroids and/or long period asteroids are more likely candidates for 
tumbling rotations - as is the case for Toutatis - since their damping timescale is 
close to the age of the Solar System. On the other hand all asteroids larger than 
~ 1 km in diameter and with a rotation period of less than about 10 hours should 
have their spin axis aligned with the largest moment of inertia axis after 100 million 
years. This damping hence explains that the vast majority of asteroids have single 
periodic lightcurves, and exhibit an invariable spin direction aligned with - or close 
to - their shortest axis of figure and also their axis of angular momentum (We do 
not mention long period precession/nutation which are not easily accessible with 
present ground-based photometry). 
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4 Figures of Equilibrium 

As seen previously, we know that (disruptive) catastrophic collisions can result in 
a rubble-pile asteroid, i.e. gravitationally bound re-accumulated aggregates with no 
internal cohesion [14]. It has been shown in [15] that asteroids could be fractured 
but still gravitationally bound without losing their integrity from the centrifrigal 
forces. Here wc arc interested in some particular class of asteroids that excludes 
the few "giants" bodies (Ceres, Pallas, Vesta,...), and the smallest ones which are 
supposed to be fully cohesive rocks. The limits in size of the rubble-pile category 
arc not well-defined but could be in the range 1km to 300km [16]. Pravcc & Harris 
[17], noting that fast rotators can hardly have negligible tensile strength, claim that 
bodies larger than « 0.15 km in size are cohesionless rubble-piles (see Fig. 4). Also 
Britt et al. [18], from an analysis of known asteroid densities, defines two categories 
of shattered objects among such rubble-pile Eisteroids. These authors distinguish the 
fractured or heavily shattered objects with porosity in the range 10-25% from the 
loosely consolidated rubble-piles with porosities in the range 30-80%. Last, a classifi- 
cation that considers the relative tensile strength together with the porosity has been 
proposed [19]. Since such a rubble-pile asteroid should be a cohesionless and gravity- 
dominated body, one can expect it to have a figure of equilibrium. Isaac Newton, 

back in 1687 in his Principia, derived the flattening (e = 1 polar radius — ^£ 

' ° ^ equatorial radius ' 

the Earth by considering it as a fluid of constant density and equalizing the weight 
of the water as due to gravitational and centrifugal acceleration in two radial canals, 
one directed toward the pole and another directed toward the equator. He could so 
explain from his theory of gravitation that the equator is not more submerged by 
the oceans than the poles. The flattening derived by Newton® e = 1/230 is that 
of the rotating Earth considered as an incompressible fluid, yet different from the 
modern value 1/298.3. Following the work of Newton other mathematicians studied 
the figures of equilibrium of a rotating mass, the reader is referred to one of the 
most comprehensive work that was made by S. Chandrasekhar [20]. A few decades 
after the work of Newton, Maclaurin (1742) in England showed that ellipsoids of 
revolution are equilibrium figures of homogeneous mass of fluid in rotation, and A. 
Clairaut (1743) in France gave a relation between the density distribution inside the 
Earth and its flattening at its surface. This was followed one century later by the 
result of Jacobi (1834), who showed that there also exist a class of tri-axial ellip- 
soids for the flgure of hydrostatic equilibrium. We will see three approaches to this 
problem considering different rlieology: incompressible fluid, compressible material 
in the linear-elastic regime, and plastic-elastic material before yield. 



4.1 Hydrostatic equilibrium 

For an incompressible fluid at rest to be in equilibrium the pressure and centrifugal 
force must balance the gravity. The equation of hydrostatic equilibrium states that, 
for each element of volume, the external body force and the boundary surface force 
must balance; so that considering the force per unit mass / and the pressure p 

® Ch. Huygens in 1690 had a very similar approach but fundamentally different 
in that he did not believe in Newton's gravitation law and his calculation yield 
£ = 1/578, i.e. as if all the mass of the Earth were concentrated in it center. 
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we have pt = Vp. When the external force represents a scalar potential field f = 
gradU and the density is constant one finds pU — p = 0, hence under gravity only 
p = |7rGp^ (-Rq ~ r^)- Similarly, stating that each volume element of the fiuid is 
at rest also results in an isotropic stress tensor Oij — —pSij where p = —ctu/S is 
the hydrostatic or mechanic pressure: i.e. all normal (and compressive) stresses are 
equal and the shearing stresses are zero. 

Introducing uniform rotation along one axis, the force still is obtained from a 
potential, so that one sees from a Vp = pVU that the external surface must also 
be an equipotential surface, with equal density and pressure. This is a necessary 
condition for figures of hydrostatic equilibrium; to be sufficient the total energy 
(gravitational, kinetic, tidal, ...) has to be minimized. For instance, considering a 
flattened sphere of eccentricity e - together with fact that total mass M and angular 
momentum J are conserved - one finds the figure of equilibrium by minimizing 
the energy E = W + T (sum of gravitational and kinetic energy) over the two 
free parameters that are e.g. the density p and the ellipticity e: dE/dp = and 
dE/de = 0. By introducing the (diagonal) potential-energy tensor: 



-L 



ly dXi 



and the inertia tensor lij = pxi Xj dx, Chandrasekhar derives the virial relation 
in tensor form: 



Uu + 111 = U22 + h2 = U33 



For a homogeneous ellipsoid of semiaxis (oi > 02 > 03), spinning along its shortest 
axis, the gravitational potential and kinetic energies are respectively: 



3 GM v-^ 2 J rr 5 J 
W = — — > Ai tti and T = „ =- 



10 R^- ^ ' 2M{ai+a. 



where Ai,A2 and A3 arc calculated in terms of Jacobi integrals involving the ellip- 
soid's ratio only (see e.g. [20]): 

/•oo 

Ai = {a2/ai){a3/ai) / ((oi/ai)^ + w) du ; i = [1..3] 

Jo 

A = ^/l + u+ ((o2/ai)2 + u) + {{az/aiY + u) (6) 

Note however that the adimensional Ai coefficients are defined differently to e.g. 
[21]; here we have Ai = 2. 

We only give a brief and general outline of the problem resolution and analysis 
of the configuration's dynamical and secular stability. The reader will find a very 
concise treatment in [20]. Here we are most interested in the figures the rotating 
fiuid can take as a the result of the hydrostatic equilibrium hypothesis. 

If there is no rotation, the equipotentials are spherical and so is the figure of 
equilibrium. A mass of homogeneous fluid with a relatively small ratio of rotational 
to gravitational potential energy T/]iy] will resemble a Maclaurin spheroid (Fig. 5), 
i.e. a axisymmetric spheroid with oi = 02 > 03 > 0. Then by putting the ellipticity: 

e = - (a3/ai)2 (7) 



we have: 
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Ai = Ai 

^1 = ^ ^sin~^ e — e\/ 1 — 
^3 = ^1 , ^ ^ - sin~^ e I 
and the rotation frequency is given by: 

By increasing for instance the rotational energy, the fluid will evolve along this 
sequence through flatter conflgurations. At a sufficiently high r/|M^| = 0.1375, and 
fP" I^-kGp) = 0.374, the equilibrium is secularly unstable, and one finds that the 
axisymmotric configuration is no longer the lowest energy state available. There is 
another sequence consisting of tri-cixial ellipsoids (ai > ai> az) with some specific 
relation ${a\la'2^ailaz) = 0, the Jacobi sequence, along which a uniformly rotat- 
ing incompressible fiuid will now evolve'^. Another sequence at still larger T/|VK| 
that bifurcates from the Jacobi sequence at some instability is one yielding binary 
structures. The latter could explain the formation from a catastrophic collision of a 
binary asteroid system through rotational fission. We will focus on the result of Ja- 
cobi that is of particular interest here, since we know from asteroids lightcurves that 
their shapes are not ellipsoids of revolution but better approached by tri-axial ellip- 
soids. So, we can suggest that the shape of rubble-pile asteroids, as re-accumulation 
of a large number of aggregates with no internal cohesion, could mimic that of 
incompressible fluids in hydrostatic equilibrium and, depending on their angular 
momentum, could be Jacobi ellipsoid. 

One possible application of such result is that knowing the shape and rotation 
period of an asteroid as derived from the lightcurve analysis, one can determine - in 
the framework of an incompressible fluid - its bulk density®. This of course assumes 
that the asteroid shape is close to an ellipsoid and that its density is constant. We 
will give two illustrative examples. Let us first remind the relation on the shape of 
a Jacobi ellipsoid: 

2 2 A\ — A2 2 . 

0-2 - 0-i 

or equivalently 

Ai - (as/ai)^ A3 = A2 - (03/02)^ A3 

Numerical results for Jacobi figures are given in 1 together with the correspond- 
ing rotational frequency sequence: 

i?' = = 2 (A2 - (aa/aa)^ ^3) 

The first example is 45 Eugenia, which body appears to have approximately the 
shape of a Jacobi ellipsoid (01/02 ~ 1.35; 02/03 ~ 1.5). Knowing that its rotation 

^ The case of viscous or non-uniformly rotating fluid, that can evolve along diflferent 

ellipsoid sequences, is not discussed here. 
* Newton originally used this scheme in his Principae to derive the density of Jupiter 
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Fig. 5. Maclaurin and Jacobi sequences. Theoretical flattening of the homogeneous 
incompressible Earth, Jupiter, Saturn and Neptune {circled crosses). The observed 
flattening of the planets are lower (crosses). We also give the sequence of an com- 
pressible spheroid with a polytrope index of n = 0.88. 



Table 1. Jacobi figures of hydrostatic equilibrium. 



1.00 0.583 0.3742 

0.95 0.598 0.3738 

0.90 0.613 0.3726 

0.85 0.630 0.3703 

0.80 0.647 0.3668 

0.75 0.664 0.3620 

0.70 0.683 0.3557 



a2/a,\ 03/02 fi 

0.65 0.703 0.3475 

0.60 0.723 0.3373 

0.55 0.744 0.3248 

0.50 0.767 0.3096 

0.45 0.790 0.2913 

0.40 0.814 0.2696 

0.35 0.839 0.2443 



02/01 03/02 fi 

0.30 0.865 0.2149 

0.25 0.891 0.1813 

0.20 0.918 0.1436 

0.15 0.944 0.1027 

0.10 0.968 0.0605 

0.05 0.989 0.0219 

0.00 ^1 ^0 
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period is P = 5.7 hr we find a bulk density of pb = 1.24. Further, assuming that this 

S-typc asteroid is constituted with material of density pg ~ 2.3 — 3, we find a macro- 
porosity of ~ 45 — 60 %. Such porosity also seems in good agreement with our rubble- 
pile hypothesis, and would roughly correspond to a random packing of aggregates. 
Moreover Eugenia is known to possess a satellite (would it be the outcome of a 
catastrophic collision, it would also enforce the rubble-pile hypothesis), that orbits 
its primary in ~ 4.7 days, so that one can independently estimate the bulk density 

+ 0.6 

Pb = M/V ~ 1.2-0-2 [22]. This is in good agreement with the value obtained from 
the hydrostatic fiuid model. This suggests that Eugenia is a Jacobi ellipsoid, and its 
overall shape on the macroscopic scale follows equipotential surface. Moreover 
Eugeiua could be a homogeneous body with a constant density profile. The second 
example is 63 Ausonia which shape is well approached by an ellipsoid but neither a 
Maclaurin (oblate) spheroid nor a Jacobi ellipsoid. However, the shortest axis being 
not very well determined by present high resolution observations, let's assume that 
the shape of Ausonia is close to that of a Jacobi ellipsoid. Then the large fiattening 
of this body (ai/ag ~ 2.2) would provide a density of pb = 0.6. Further, assuming 
this C-type asteroid is constituted of material with density pg ~ 2 — 2.5, we find 
a macro-porosity of ~ 70 — 76 %, which value seems rather unrealistic! On the 
other hand, the observed lightcurves amplitude of Ausonia could be obtained with a 
binary asteroid where each component, in hydrostatic ecjuilibrium [23], would have 
a somewhat higher density [24]. We shall however see in Section 5.3 that Ausonia 
actually is a single body with a shape close to a prolate spheroid. Last, if one plots 
the observed asteroids shapes against the Maclaurin and Jacobi shapes (see Fig. 6), 
one clearly sees that asteroids shapes generally departs from equipotential surfaces. 

If Eugenia cam be considered as a Jacobi ellipsoid, it seems nevertheless that, 
as seen before, asteroids in general and rubble piles in particular do not follow such 
figures of equilibrium for inviscid and incompressible fiuids. Interesting, instead of 
the Jacobi ellipsoid shapes, the prolate spheroid (ai > 02 = as) appear to be a more 
common shape among asteroids. The previous analysis can be completed, in the 
case of stars as well as in the case of solar system bodies, by considering - inviscid 
- compressible fluids. In this case the density and pressure are no more constant 
through the body, and are related by some general law f{p,p,T) — 0. For instance 
Laplace considered the relation dp/dp = hp, and Roche dp/ dp = hp + h' p^, h and h' 
being constants. In a similar way one can also consider bodies of incompressible but 
non homogeneous fluids. In addition to the numerical experiments simulating the 
behavior of compressible stars, let us mention the analytical work achieved in [25] for 
describing the ellipsoidal figures of equilibrium in the compressible case. The authors 
have considered polytropes of index n,p = Kp^ with F = 1 + 1/ n. Briefly, such fluid 
compressibility will change the shape of the sequence as can be seen in Fig. 5, but 
not the overall shape of the flgure of equilibrium. In other words, Maclaurin spheroid 
and Jacobi ellipsoids are still figures of equilibrium but the rotational frequency for 
a given shape is different to that of the incompressible case. Moreover, as shown 
in Fig. 5, for a given shape the rotational frequency is larger when considering this 
particular density distribution; thus the bulk density for a given shape is smaller 
than it is in the incompressible case. From that it also appears that the knowledge 
of the rotation period and geometric flattening (or shape) is not sufficient to obtain 
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Fig. 6. Observed asteroids axis ratio against Maclaurin spheroids and Jacobi el- 
lipsoids. Are overplotted the data for Eugenia on the right {filled squares), and on 
the left the data for the prolate spheroid Ausonia as well as the hypothetical Jacobi 
shape Ausonia. 

information on the density distribution inside the body. Also such compressibility 
would not provide a more realistic bulk density for our Ausonia example. 

4.2 Elastostatic equilibrium and elastic-plastic theories 

Another approach of interest for solid bodies and that extends the approach of fluids 
is to consider elastic bodies^. We now consider a rheology where there is a linear 
relation between stress aij and strain etj, and the law of constraint is given by 
Hooke's law : 

aij = A ekk 5ij + 2/i e^j 

where (A, ^) are the two Lame coefHcients of elasticity. The strain is derived from 
the deformation field Ui(xi) in the frame (xi,X2,X3) by: 

1 / dui duj \ 
" 2 \dF, ^Ih^) 

and the srain-stress relation can also be written as (see e.g. [26]): 
^ And, say, thermoreologically simple bodies 
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Eeij = + (Tij - vJ^k'^kkSij (10) 

where -B > is Young's modulus, and —1 < v < 1/2 the Poisson^" ratio, coefficients 
that only depend on the considered material. For instance, in the case of a simple 
compression an < 0, 1722 = ctss = along axis xi, one has en = l/Ean, and 
the corresponding deformations rate along the perpendicular directions are £22 = 
633 = —vjEow — — i^eii. The couple of parameters {E,v) is uniquely connected to 
the two Lame coefficients. These coefficients also express the mechanic dissipation 
inside the material. Considering by continuity the limiting case where v —> 1/2, 
the shearing-stiffness modulus is n = and one should find the results for the 
incompressible and inviscid fluids case. The static equilibrium is obtained when the 
sum of internal and external body forces vanishes div 51 + f = 0, whore 51 — [(Tij] 
is the stress tensor. Considering the gravitational and centrifugal potentials from 
which the forces are derived (f = pgrad(U)), one obtains: 

As seen above in Sect. 3 we can assume that the object is spinning along its shortest 
axis with constant angular rate i? and write the centrifugal potential 

Uc = {xi + xi)/2 

The gravitational potential is given by: 

UG = -KGpa\{Uo-^Ai xl) 

i 

where the coefficients have been defined in (6), and (7o = du. Solution of 

equilibrium can next be obtained by minimization of the elastic energy. The problem 
is analytically tractable because the loads are linear in the spatial coordinates and 
because of the symmetries in the considered figure. Por instance this approach has 
been applied to the non-rigid spheroidal Earth [28]. The analytical treatment in the 
case of tri-axial ellipsoids is more cumbersome and has been treated in e.g. [29, 30, 
31]. Considering an homogeneous, isotropic, linear-elastic, and slightly compressible 
material with v = 0.45; 0.499, Washabaug & Scheeres [31] have shown that, at low 
angular momentum, ellipsoidal figures with compressive stresses at the surface exist 
(i.e. the presence of tensile strength is not needed) but they generally lie only in 
the vicinity the minimal energy state. They also showed that the elastic energy 
minima are weaker for compressible material, but nevertheless they occur at high 
angular momenta and for very elongated shapes which shapes are not observed 
among asteroids. 

Eventually, wo will discuss a more general approach that does not depend on the 
actual stress-strain behavior or possible residual stresses, but we will consider the 
limit stresses of an elastic-plastic body. Starting from (11) - and considering only 
ellipsoids - the equilibrium equations and boundary conditions leaves three degree 
of freedom in the general solution. Thus one can consider that the material is fluid 
(imposing that the shear stresses are vanishing and all normal stresses are equal). 

In practice we have i/ > although negative Poisson's ratio have been witnessed 

in some foam [27], see URL http://sllver.iieep.wlsc.edu/~lcLkes/Poisson.html 
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or one can consider linear-elastic isotropic material by introducing an additional 
relation between strain and stresses from Hooko's law. These arc the two options we 
have considered so far. Considering fluids or elastic deformations to model asteroids 
would not allow for large fails, boulders or craters on the surface of these small rocky 
bodies as can be witnessed for instance on the surface of the planetary satellite 
Phobos or asteroid (253) Mathilde. It has been suggested that cohesionless bodies 
could maintain shapes significantly different from figures of hydrostatic equilibrium, 
so long as i) one assumes that the rotating mass is not a fluid but behaves as a 
granular soil with non negligible solid-to-solid friction, and ii) internal stresses are 
not high enough to crush individual particles [32, 33]. In this case the material 
can sustain non negligible tangential force before plastic flow. The Mohr- Coulomb 
criteria is generally used for deriving the maximal stress strain before yield of a 
given soil. The stress tensor II = [aij] being symmetric, it is diagonalizable; given 
the principal stresses ai, a2, and neglecting the cohesive strength, the Mohr- 
Coulomb criteria states: 

, , ^ 0-1/0-3 - 1 
tan > — 

2-yiTl/fT3 

that is, the maximal tangential stress is limited by only the most and least (com- 
pressive) normal stresses sigmaiandaa. In such case we no more have an equality 
that provides the unique figure of equilibrium, but instead a range of possible con- 
figurations. Considering again tri-axial ellipsoids of constant density, it appeaxs that 
the ratio sigmai/as does not depend on the spatial coordinates, and one can plot 
the rotational frequency versus the axis ratio for a given friction angle [33]. We see 
in Fig. ?? that the hydrostatic equilibrium is obtained for negligible friction, and for 
a given friction angle (e.g. # = 15deg) there is a zone of possible ellipsoidal figures. 
That is, for those tri-axial ellipsoids the gravitational and centrifugal load are small 
enough to be sustained by the friction and hence avoid failure. We also see that the 
rotation frequency is limited at the largest possible friction $ = 90 deg, in particular 
for a sphere one has [4]: 

V ^max = 47r/3 

but this does depend linearly on the axis ratio. Last if one plots the data from 
known asteroid, assuming typical densities in the range 1-3, it eventually appears 
that the great majority of observed shapes is consistent with such a cohesionless 
Mohr-Coulomb model and a friction angle of « 25 deg, which value seems realistic 
and typical of dry terrestrial soils (it is « 30 for sand) . 

4.3 Binary systems and the density profile 

Asteroids do have satellites too! Since the first one that was discovered by the Galileo 
probe during its fly-by with Ida in 1993, about a dozen additional systems have been 
identified by ground-based observations (radar, adaptive optics, photometry) in the 
NEO or KBO population as well as in the main belt. Note that previous attempts 
from the Hubble Space Telescope imaging instrument were unsuccessful [34]. The 
presence of a satellite (generally 10-20 times smaller than the primary) is of high 
value to determine the mass of the primary from Kepler's third law and eventually 
its bulk density. Moreover, since these objects are rather fiattened or elongated, 
the dynamical perturbations due to the non-spherical gravitational potential are ex- 
pected to be relatively large. As seen previously the (secular) perturbations depends 
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on the potential, or mass distribution [35]. If the spin rate of the primary is much 
larger than the orbital rate of the satellite (i7rot 2> n), the secular effect of the C22 
is negligible, and the major perturbation arises from the J2 coefficient: 

, ^ h+h- 2/3 
J2 = = ^^^^^ 

Further, the secular terms are obtained from Lagrange equations: 
a = e = i = 

(l-5cos^i) C20 + - 



4a2(l-e2)2 

o „ „2 



^ - L2 (1-12)3/2 (1 - 3 cos^ C20 + ... (12) 

We know since Clairaut and Radau that the knowledge of the dynamical and 
geometric flattening brings insight on the density distribution inside the body. Thus 
observing the secular perturbation on the satellite's orbit over successive months 
provides the dynamical flattening while observing photometric variation over suc- 
cessive apparitions provides the geometric flattening. By comparing both one can 
at least test the hypothesis of constant bulk-density inside the (primary) aster- 
oid. This has been applied to the orbit of Kalliope's companion [36]. Assuming a 
homogeneous primary with constant density distribution the observed geometrical 

flattening provides a dynamical J2 = jq B.nd in turn a precession rate of 

w ~ 0.3 deg/day. This value is in severe conflict with the observed one ~ 0.7 deg/day. 
As shown in Fig. 8 there are three possible explanations that maybe all concur to- 
gether to this discrepancy: the size of the body, the geometrical flattening, and the 
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density distribution. Asteroids diameters are essentially given by an indirect method 
from observations with the IRAS satellite, and can be in error by 10% or more [37]. 
The shape being derived from disk-integrated photometric data and not high res- 
olution imaging is not known with the best accuracy neither. Last, including an 
empirical but simple density distribution of the form p{r,9, 0) = po r"; a & M, i.e. 
which increases as we progress toward the surface (r = 1), one can write the zonal 
harmonic as a function of the one for a = 0, i.e. at constant density 7° : 

5 a + 3 
3 a + 5 

Hence for a given shape, the pericenter precession is increased by concentrating the 
mass at the outer surface of the body. Assuming that the large precession is due 
to the non homogeneous mass distribution alone, one finds a density at Kalliope's 
surface of w 7, which is marginally acceptable. Let's now consider that Kalliope is a 
size-sorted rubble-pile, where the larger and more irregular rocks are in the central 
part and the smaller material is kept by friction at the outer layers [38]. Such mass 
distribution could correspond to a body of homogeneous but size-sorted material 
with larger relative voids (or porosity) in the central part, and more densely packed 
material toward the surface. 



5 The determination of shape and spin parameters by 
Hubble Space Telescope 

5.1 The FGS interferometer 

The Hubble Space Telescope (HST) is a complex of instruments, built around a 
2.4-meters telescope, orbiting the Earth. For several reasons (pressure on observing 
time requests, time constraints due to the orbit, complexity of the instruments, etc.) 
it is sometimes very difficult to apply for successfully and to use. 

However, its great advantage - being outside the Earth atmosphere - allows 
to obtain otherwise difficult measurements and observations. While its imaging ca- 
pabilities of deep-sky objects and planetary surfaces/atmospheres are well known, 
they are not sufficient to perform accurate measurements of asteroid shapes and 
sizes. In fact, the highest resolution is reached by the Planetary Camera, having a 
plate scale of 46 mas/pixel^ ^. This value is of the same order than the apparent size 
of several, interesting main belt objects at opposition, and allows some resolution 
to be achieved on a very restricted set of the largest bodies only. 

A much more sentive instrument is the Fine Guidance Sensor (FGS) , an interfer- 
ometer normally used to allow careful pointing and guiding of the HST while imaging 
is performed by the main CCD cameras. Three FGS instruments are mounted close 
to the focal plane, and each of them works by producing interference between the 
two beams coming from the defocalised semi-pupils of the telescope. Inside the FGS 
the beam is divided into two parts, associated to two perpendicular axis (in the fol- 
lowing we refer to them as the FGS-X and FGS-Y axis). Each beam enters a Koester 

In the following we express all apparent sizes in milli-arcseconds (mas) = 10~* 
seconds of arc 
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Fig. 8. Secular periastron advance CJ20 for the orbit of 22 Kalliope's companion. 



prism. Inside this device the self-interference occurs. The resulting two beams ex- 
iting a single Koester are then collected by two photometers, measuring their flux. 
The important feature to note here, is that the difference between the fluxes is a 
function of the inclination of the incoming wavefront^^. 

The FGS is able to "scan" the focal plane in straight line, along the FGS-X and 
FGS-Y axis, in steps of 1-2 mas, and over a few arcseconds. Each step corresponds 
to a different inclination of the wavefront and thus to a different response of the 
photometers. A response function, called "S-curve", is reconstructed from the flux 
difference normalized to the total flux. 

It is clear, then, that the response for an extended source depends upon the 
light distribution on the focal plane. The FGS sensitivity being optimal for spatial 
frequencies < 200 mas, it is normally employed to measure the diameter of large 
star disks or the separation of close doubles. An example of the simulated response 
curve for an extended, uniform disk, or for a double disk, is given in figure (Fig. 

More details can be learned by browsing the HST/FGS on-line handbook at URL: 
http : / /www . stsci . edu/ instruments/ f gs/. 
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Fig. 9. Upper panels: the simulated response curve, along the FGS-X and FGS- 
Y axis, for a uniform, circular source of different sizes (indicated in arcsec). The 

"template" curve correspond to the observation of a point-hke star. Lower panel: 
the simulated response curve for a double source composed by two equal, tangent 
disks of different diameters. 
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9). As it can be seen, the largest is the source, the smallest is the amplitude of the 

S-curvc. A double object has the tendency to double the curve peaks. In order to 
better understand the results presented in the following, it must be noted that each 
scan by the FGS produces two S-curves, one for each axis. As a consequence, the 
extracted informations concern the size of the studied object as projected along the 
two FGS axis. 

In recent years, the FGS has been used for the first time to measure rapidly 
rotating and moving objects: the asteroids. Six objects having interesting spin and 
lightcurve properties have been chosen, each observed for an HST "orbit", i.e. for 
a total duration of about 40 minutes (corresponding to 30-40 FGS scans). During 
this time lapse two movements affect the observation: the spin of the object and its 
proper motion. The first one is highly desirable, since it can help to constrain the 
asteroid shape by studying the variations of the FGS-X and FGS-Y projections. The 
second one, on the other hand, has been minimized by choosing epochs close to the 
stationary points in the asteroids apparent motion, and by applying a correction (a 
posteriori) to compensate for the parallax due to HST displacement. 

In summary, each set of scans for each asteroid contains informations on its 
shape and size, as projected over FGS-X and FGS-Y axis, over a limited fraction of 
the rotation period - roughly corresponding, for the given set of objects, to about 
the 10% (about 30-40 degrees of rotation). It can thus be supposed that a fitting of a 
model to the set of S-curves can lead to the reconstruction of sizes and shapes. The 
degree of accuracy (i.e. the number of parameters to be determined) will depend 
upon the quantity of observations and their signal/noise ratio. 

5.2 From data to modeling 

In general, given a set of FGS-X and FGS-Y S-curve couples at each epoch, the 
fitting procedure begins by determining shape, orientation and size of the on-sky 
projected ellipse that separately accommodates best to each curve. 

The problem, then, is to identity which is the solid, three-dimensional body that 
is capable of producing the observed, projected ellipses. If a traditional equilibrium 
figure is searched for, an ellipsoid with a given orientation of the spin axis will be 
the figure of choice. 

It must be noted that traditional photometry from Earth-based telescopes nor- 
mally records brightness variations associated to the object shape, directly yielding 
its rotation period. From lightcurves taken at different epochs, a first estimation 
of the spin axis direction can be obtained. In general, however, some symmetries 
in the problem does not allow to discriminate between degenerate solutions. The 
result of photometric pole determination is thus expressed by two possible spin axis 
directions, each yielding, for a given epoch, the same object area projected on the 
sky, i.e. the same brightness. However, an instrument capable of directly detecting 
the orientation of the shape and its variation in time, can immediately eliminate the 
ambiguity and help to discard one of the two solutions. This is the first result that 
have been derived by HST observations. 

Thus, having selected the good pole solution between the pair available, the 
model is fitted to the S-curve is that of an ellipsoid of uniform brightness. The 
shape is described by the three ellipsoid semi-axis a, b and c. Since the uncertainty 
in pole coordinates can reach several degrees, a trial and error adjustment, reducing 



20 



D. Hestroffer and P. Tanga 



the 0-C, is performed. In fact, we realized that additionally solving the fit for pole 
coordinates seems not to add significant improvements given the sot of data currently 
available. Thus, pole coordinates do not enter the core of the fit process. 

After having determined the best-fitting ellipsoid, we will have a set of values 
corresponding to semi-axis sizes, spin pole coordinates and rotational phase. In turn, 
this allows to re-compute for each epoch the ellipse projected on the sky. Another 
iteration can thus be performed repeating the whole process from the beginning. If an 
ellipsoidal shape solution exists, the parameters rapidly converges and the residuals 
collapse. Significant residuals are, on the contrary, the indication of some detectable 
departure from the simple assumption on the shape, and some more complex models 
deserve to be taken into account. 

5.3 Some significant examples 

In the following we illustrate some results that we think to be particularly significant. 
The accuracy of the results, their limitations and the practical difficulties should be 
apparent. A complete review of results is published in [39, 40, 41], while the following 
table reviews the main parameters derived. 

Table 2. A summary of FGS size measurements, obtained by considering ellipsoidal 
models. (216) Kleopatra, having a complex shape, is not shown here. The last column 
gives the ratio of the axis; the parentheses indicate that either 6 or c are not well 
constrained. 



Name 


a, b, c [km] 


a/b 


a/c 


(15) Eunomia 


181, 103, 


102 


(1.76) 


1.78 


(43) Ariadne 


45, 26, 


26 


1.71 


(1.71) 


(44) Nysa 


59, 35, 


35 


(1.72) 


1.72 


(63) Ausonia 


75, 33, 


33 


2.28 


(2.28) 


(624) Hektor 


62, 28, 


28 


2.21 


2.21 



63 Ausonia 

63 Ausonia was one of the brightest objects in the set, and the first to be observed. 
The signal-to-noise ratio is rather good, so the fitting process operates in ideal 

conditions. Figure 10 shows the S-curves of (63) Ausonia for a selected epoch. The 
residuals of the fit with a three axis ellipsoid are shown in the inset and are very 
small in this case, showing that the identified 3D ellipsoid is completely consistent 
with the available FGS data, as discussed above. 

Due to its orientation on the sky, the asteroid, while rotating, exhibits some 
interesting variations in the length of the projected a and h axis, while c (coinciding 
with the rotation axis) does not move and its infiuence on the FGS measurements 
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Fig. 10. Tlie Ausonia S-curves at tlie beginning of the observation, fitted by using 
the tri-axial ellipsoid given in Table 2. The residuals of the fit given in the inset are 
very small. 



is minor. For this reason, the value of the c axis is affected by a high uncertainty 
(around 5 mas or more) while a and b are constrained to a 1 mas level. These figures 
can be considered to be typical for this kind of model-dependent fit. 



216 Kleopatra 

This was most irregular body observed. A little time before the HST observations, it 
was observed by radar [42] . The reconstructed shape seemed to hint to a bi-lobated 
object, very elongated and irregular. 

The HST/FGS signal suggests an elongated shape, well approximated by two - 
not detached - ellipsoids, whose signature is well visible in the S-curve. The overall 
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shape appears to be more elongated and flattened in comparison to radar data. 
Details are given in [40] and [43]. 
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Fig. 11. The variation of the sizes projected on FGS-X and FGS-Y for the main 
component of the double-lobed asteroid (216) Kleopatra, during the observation. 
The shape of the object and its orientation relatively to the FGS axis are represented 
by the two ellipsoids below the curves. 
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Fig. 12. Tlie indicative sliape of Hektor obtained by a best-fit of tfie FGS observa- 
tions. Tlie noisy data set does not allow to draw definitive conclusions. 

624 Hektor 

This Trojan asteroids is the faintest asteroid observed by HST/FGS during the 
program (V=i5.0 at the epoch of the observation), being at an average distance 
corresponding to the semi-major axis of the orbit of Jupiter. Due to its faintness, 
the S/N ratio is very small and it probably represents the lowest possible for this 
kind of studies by the FGS. The response functions is best fitted by a very elongated 
shape, but due to the low S/N ratio it is not possible to clarify if (624) Ifektor can 
really be considered to be a contact binary as supposed in the past by Weiden- 
schilling [44]. Figure 12 shows the best-fit shape obtained by HST observations. 



15 Eunomia 

Together with the previous object, (15) Eunomia is another difhcult target. (15) 
Eunomia is the object having the largest apparent size, as shown by the small 
amplitude of the S-curve. The shape of the signal is consistent with a shape more 
complex than a simple ellipsoid, as shown in Fig. 13. Unfortunately, in order to 
better constrain its shape, more observations and a widest coverage of its rotation 
would be necessary. 

6 Conclusions 

The observations by HST/FGS are a sensitive and powerful method to determine 
size and shape of asteroids. However, some limitations should be kept in mind. 

First of all, the observing time is very limited and hard to obtain. For this reason, 
it is very difficult to sample the whole rotation curve of an object. Unfortunately, 
this is the only way to constrain complex shapes. 
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Fig. 13. Lower panel: the suggested shape of (15) Eunomia is consistent with the 
strongly asymmetric S-curve obtained for this object. Other possible solutions re- 
quire the use of other constraints, such as those coming from photometry. 



Furthermore, the shape solution that is found depends upon some a priori 
choices, such as the basic shape model, its scattering properties, the absence of 
albedo markings, etc. This details, while not changing the qualitative interpretation 
of the data, are probably important in order to define the ultimate precision of the 
observations. 

Finally, the FGS can work on a limited sample of extended and bright asteroids. 
However, the same basic approach, coupled with more sophisticated data inversion 
techniques, can be applied to other optical interferometers, based on the Earth 
surface. We can thus hope that, in the future, the sample of asteroids for which 
size, shape and orientations are known will give us a more complete view of their 
physical properties. Last, combining photometric data to high resolution data and 
astrometric positions of asteroid satellite should bring insight of the primary asteroid 
interior and possible coUisional history. 
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